SLAC-PUB-6576 
July 1994 
(T/E) 



Polarization as a Probe to the Production Mechanisms of Charmonium in vrA^ 

Collisions^ 

W.-K. Tangg and S. J. Brodsky^ 

Stanford Linear Accelerator Center 

'nJ" ■ Stanford University, Stanford, CA 94309 

G\ ■ 

a^ : 

I , M. Vanttinen and P. Hoyer 
bJOl Department of Physics 

^ I University of Helsinki, Finland 

^ ; 

(N . 

Abstract 

T— I ■ 

I Measurements of the polarization of J/V' produced in pion- nucleus collisions are in 

' disagreement with leading twist QCD prediction where J/ip is observed to have negli- 

. gible polarization whereas theory predicts substantial polarization. We argue that this 

OO I discrepancy cannot be due to poorly known structure functions nor the relative produc- 

^ ' tion rates of J ftp and XJ- The disagreement between theory and experiment suggests 

■ important higher twist corrections, as has earlier been surmised from the anomalous 

I non-factorized nuclear ^-dependence of the J/i/j cross section. 

1 Introduction 

; 

>• ■ One of the most sensitive tests of the QCD mechanisms for the production of heavy quarko- 

^ ! nium is the polarization of the J/ip in hadron collisions. In fact, there are serious disagree- 

ments between leading twist QCD prediction and experimental data |^, ||, |], ||, 0] on the 
production cross section of 'direct' J/ip and xi- We would like to advocate that polariza- 
tion of J/ijj provides strong constraints on the production mechanisms of J/ip and thus can 
pinpoint the origin of these disagreements. 

In this paper we will present some preliminary results on the theoretical calculation of 
the polarization of J/ip in vrA^ collisions. The completed analysis will be published in a 
later paper [|^. We found that the polarization of J/ip provides important constraints on the 
nature of the production mechanisms and urge that polarization measurement of J/ ip should 
be included in the design of future charm production experiment. 

The paper is organized as follow. In section 2, we show that from the experimental data 
on the production cross sections and leptonic decay widths of direct J/ip and ip', the long 
distance physics of formation of bound states of cc can be separated from the short distance 
physics of production of the cc pair. Thus, the perturbative analysis is under control in 
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calculating J/ip production even though the mass of charm quark is not much larger than 
^QCD- Once the validity of perturbative method is established, we calculate the production 
cross sections of direct J/ip., Xi ^"^^ X2 in ttA^ collisions in PQCD. These results are presented 
in section 3 and discrepancies are observed. We show that, in comparison with the recent 
E705 and E672 data p, [10|, the predicted ratio of direct J/ip production compared to the X2 
production is too low by a factor of about 3. In addition the production ratio of production 
cross sections of xi to X2 is too low by a factor of 10 compared to data. A similar conclusion 
has been reached in |]rT| . The polarization data of J/ ip [0, |13|, Q allows us to make further 
conclusion of the origin of the disagreements. In section 4, we find that even if the relative 
production rates of the J/ip, Xi and X2 are adjusted (using i^-factors) to agree with the 
data, the J/ip polarization data is still not reproduced. Therefore, the discrepancies do not 
arise from an incorrect relative normalization of the various channels and new production 
mechanisms are needed. We will present our conclusion in the last section. 



2 Can direct J/ip production be calculated in PQCD? 

In leading twist QCD, the production of the J/tp at low transverse momentum occurs both 
'directly' from the gluon fusion subprocess gg —>■ J/ip + g [Fig. la] and indirectly via the 
production and decay of Xi and X2 states. These states have sizable decay branching fractions 
Xi,2 ^ J/ip + J of 27% and 13%, respectively. 



Figure 1: Fig. la shows direct J/ip production through gg scattering. The formation of bound 
state is described by the wavefunction \Ef}/^(0) at the origin. Fig. lb shows leptonic decay 
of J/ijj into e~^e~ pair. The probability of finding the cc pair is given by the wavefunction 
*J/^(0). 

In this model, we assume that the non-perturbative physics, which is described by the 
wave function at the origin in cases of production of J/tp and ip', is separable from the 
perturbative hard subprocess, i.e., factorization holds. As the wave function at the origin can 
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be related to the leptonic decay amplitude [Fig. lb] , the ratio of ip' to direct J/ ip production 
can be expressed in terms of the ratio of their leptonic decay width. More precisely, taking 
into account of the phase space factor, 

a(ij') TO;' e+e") Mj, 

> ^ > ■'1^ ~ 0.24 ±0.03 1 



ad^r{J|'lp) r(J/^ ^ e+e-) M^, 

where adiriJ/ip) is the cross section for direct production of the J/'ip- The ratio (|l]) should 
hold for all beams and targets, independently of the size of the higher twist corrections in 
producing the point-like cc state. The energy should be large enough for the bound state 
to form outside the target. The available data is indeed compatible with ([^). In particular, 
the E705 value [§] is 0.24. In Table 1, the ratio of ip' to direct J/ip production with different 
projectiles is presented. They are all consistent with the value 0.24. 





aii^') [nb] 


CTdiriJ/tp) 




71+ 


22 ±5 


97± 14 


0.23 ±0.07 


7T^ 


25 ±4 


102 ± 14 


0.25 ±0.05 


P 


20 ±3 


89 ± 12 


0.23 ±0.05 



Table 1: Production cross sections for tp', direct J/tp and their ratio in tt+A^, vr and pN 
collisions. The data is from Ref. [§. 

The anomalous nuclear target A-dependence observed for the J/ip is also seen for the 
ip' |TB|, so that the ratio (|I|) is indeed independent of A. Therefore, at high energies, the 



quarkonium bound state forms long after the production of the cc pair and the formation 
process is well described by the non-relativistic wavefunction at the origin. 



3 Production rates of ip and xj states at leading twist 

In leading twist and to leading order in a^, J/ip production can be computed from the 
convolution of hard subprocess cross section gg — > J/ipg, gg — * Xj? cic, with the parton 
distribution functions in the beam and target. Higher order corrections in as, and rela- 
tivistic corrections to the charmonium bound states, are unlikely to change our qualitative 
conclusions at moderate xp- Contributions from direct J/tp production, as well as from 
indirect production via Xi ^ind X2 decays, will be included. Due to the small branching 
fraction xo J/ip + 1 of 0.7%, the contribution from xo to J/ip production is expected 
(and observed) to be negligible. Decays from the radially excited 2^5"! state, ip' J/ip + X, 
contribute to the total J/ip rate at the few per cent level and will be ignored here. 
The 7tN —^X2 + X production cross section to lowest order is 

cr(7riV ^ X2 + X] xf>0) = / Fg/^{xi)Fg/N{r/xi)ao{gg -> X2) (2) 

J y/r Xi 

where r = M^^/s and the quantity aQ{gg X2) = 167r^a^|i?p(0)p/M^2 |T^. We restrict 
the X2 momentum range to the forward CM hemisphere {xp > 0) in accordance with the 
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available data, and use the structure functions of Ref. WB, O evaluated at = M^^ . We 
also take the renormalization scale to be = M"^, . 

The direct vrA^ — * J /ip + X cross section is similarly given by 

a{TcN ^ J/ijj + X; xp > 0) = / dxi / dx2 / diFg/T,{xi)Fg/j^{x2) 

Agg^J/i, + g) (3) 
at 

where t is the invariant momentum transfer in the subprocess, and 

/ X2M'],^ - XiS \ 

t^^ = max '-^^ , Ml^ -s\. 4 

Eq. (|^) also applies to the vrA^ — 5> Xi + X reaction, in which case a sum over the relevant 
subprocesses gg Xi9j 91 Xili 91 ~^ X1I and qq Xi9 is necessary. The differential 
cross sections dcr/dt for all subprocesses are given in IT^, O. 



In Table 2 we compare the X2 production cross section, and the relative rates of direct 
J/ip and xi production, with the data of E705 and WAll on n'N collisions at Eiab = 300 
GeV and 185 GeV M. 





(^iX2) [nb] 


crdiriJ/i')/criX2) 


^(Xl)/^(X2) 


Experiment 


188 ±30 ±21 


0.54 ±0.11 ±0.10 


0.70 ±0.15 ±0.12 


Theory 


72 


0.19 


0.069 



Table 2: Production cross sections for xi, X2 and directly produced J/ip in tt~N collisions. 
The data from Ref. H include measurements at 185 and 300 GeV. The theoretical 
calculation is at 300 GeV. 

The X2 production rate in QCD agrees with the data within a 'i^'-factor' of order 2 to 
3. This is within the theoretical uncertainties arising from the J/tp and x wavef unctions, 
higher order corrections, structure functions, and the renormalization scale. A similar factor 
is found between the lowest-order QCD calculation and the data on lepton pair production 



[pOj , pi|] . On the other hand. Table 2 shows a considerable discrepancy between the calculated 
and measured relative production rates of direct J/ip and xi; compared to X2 production. 
A priori we would expect the i^-factors to be roughly similar for all three processes. We 
conclude that leading twist QCD appears to be in conflict with the data on direct J/ip and 
Xi production. Although in Table 2 we have only compared our calculation with the E705 
and WAll 7r~N data, this comparison is representative of the overall situation (for a recent 
comprehensive review see ||TT[|). 



4 Polarization of the J/ip 

The polarization of the J/ is determined by the angular distribution of its decay muons in 
the J/ip rest frame. By rotational symmetry and parity, the angular distribution of massless 
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muons, integrated over the azimuthal angle, has the form 



da 



dcosO 



oc 1 + A cos^ 9 



(5) 



where we take 9 to be the angle between the /i"*" and the projectile direction {i.e., we use 
the Gottfried-Jackson frame). The parameter A can be calculated from the cc production 
amplitude and the electric dipole approximation of radiative x decays. 

The electric dipole approximation of the radiative decay Xj ^ 4'1 is exact in the heavy 
quark limit; i.e., when terms of 0{E^/mc) are neglected. As a consequence, the heavy quark 
spins are conserved in the decay, while the orbital angular momentum changes. 

The lowest order subprocess g{fJ^i)g{fJ^2) ^ cc ^ X2{Jz) only produces X2 with = ±2 
states asumming that the transverse momenta of the incoming gluons are neglected. In 
the Jz = ±2 polarization state the spin and orbital angular momenta of its constituent 
charm quarks are aligned, = = ±1. Since Sz is conserved in the radiative decay 
X2 J/i^ + 7, it follows that JziJ/i^) = Sz = ±1 (L = for the J/^p). Thus the J/ip's 
produced via X2 decay are transversely polarized, i.e., A = 1 in (||). This result is exact 
if both the photon recoil and the intrinsic transverse momenta of the incoming partons are 
neglected. Smearing of the beam parton's transverse momentum distribution by a Gaussian 
function exp [— (fc^/SOO MeV)^] would bring A down to A ~ 0.85. 

From the gg ^ J/ip + g amplitude we find for direct J/ ip production, ttN J/ip + X 
fi+(x- +X, 



da 



dxidx2 

BnndxFdcos9 6477 7 (xi + X2)s 



Fg/n{Xl)Fg/N{x2) 



011 + 000 + {011 - 0oo) cos^ 9 



(6) 



where is the J/tp A*'''/^" branching fraction, xp = '^P^/Vs is the longitudinal- 
momentum fraction of the J/ip, and 9 is the muon decay angle of Eq. (|]). The f?ii,f?oo 
are the density matrix elements and can be found in 

For the ttN —>■ xi + ^ ~^ J/ip + •j + X — >• + 7 + X production process we get 

similarly 

1 da 3 ^ , > ^^ f dxidx2 
Br(xi VV^ij 



B^^ dxpd cos 9 



12871 



X 



^00 + ^011 + (000 



{Xi + X2)s 

0ii) cos^ 9 



Fii^{Xi)Fj/N{x2) 



(7) 



where the density matrix elements for ij = gg, gq gq and qq scattering are again given in 

In Fig. 2a we show the predicted value of the parameter A of Eq. (^) in the GJ-frame 
as a function of xp, separately for the direct J/ip and the xi,2 — ^ J/ip + 'y processes. Direct 
J/ip production gives A ~ 0.25, whereas the production via xi results in A ~ —0.15. 

The A(xir)-distribution obtained when both the direct and indirect J/ip production pro- 
cesses are taken into account is shown in Fig. 2b and is compared with the Chicago-Iowa- 
Princeton [0 and E537 data for 252 GeV irW collisions and 150 GeV 7t~W collisions 
respectively. Our QCD calculation gives A ~ 0.5 for xp ^ 0.6, significantly different from 
the measured value A ~ 0. 
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Figure 2: CIP (•) and E537 (o) data compared with theoretical prediction. Fig. 2a shows 
the parameter A from different contributions: direct J /ip, Xi,2 J/ip + 'J processes. Sohd 
curves shows the results with the intrinsic transverse momentum of the incoming partons 
neglected while the dashed curves have the beam parton's transverse momentum modeled 
by a Guassian function exp[— (/c^/SOOMeV)^]. Fig. 2b takes into account both the direct and 
indirect J/ip production: without K factors correction (solid curve), and with K factors 
correction (dashed curve). 

The discrepancies between the calculated and measured values of A is one further indica- 
tion that the standard leading twist processes considered here are not adequate for explaining 
charmonium production. The J/ip polarization is particularly sensitive to the production 
mechanisms and allows us to make further conclusions on the origin of the disagreements, 
including the above discrepancies in the relative production cross sections of J/ip, Xi ^i-nd X2- 
If these discrepancies arise from an incorrect relative normalization of the various subprocess 
contributions {e.g., due to higher order effects), then we would expect the J/ip polarization 
to agree with data when the relative rates of the subprocesses are adjusted according to the 
measured cross sections of direct J/ip, Xi and X2 production]^. The dashed curve in Fig. 2b 
shows the effect of multiplying the partial J/ip cross sections with the required i^'-factors. 
The A parameter is still predicted incorrectly over most of the xp range. 

A similar conclusion is reached (within somewhat larger experimental errors) if we com- 
pare our calculated value for the polarization of direct J/ip production, shown in Fig. 2a, 
with the measured value of A for ip' production. In analogy to Eq. (|I|), the 7/)' polarization 
data should agree with the polarization of directly produced J/ip^s, regardless of the produc- 

■^In the case of Drell-Yan virtual photon production, it is known that higher-order corrections do not 
change the 7* polarization significantly [ p2| , which makes it plausible to represent these corrections by a 
simple multiplicative factor, which does not affect the polarization of the photon. 
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tion mechanism. Based on the angular distribution of the muons from ip' fi^^^ decays 
in 253 GeV Tr~W coUisions, Ref. quotes A^/ = 0.02 ± 0.14 for xp > 0.25, appreciably 
smaller than our QCD values for direct J/^'s in Fig. 2a. 



5 Discussion 

We have seen that the J/ip and xi hadropro duct ion cross sections in leading twist QCD are at 
considerable variance with the data, whereas the X2 cross section agrees with measurements 
within a reasonable i^'-factor of 2 to 3. On the other hand, the inclusive decays of the 
charmonium states based on the minimal perturbative final states {gg and qqg) have been 
studied in detail using perturbation theory ||2^, ^ and appear to work fairly well. It is 
therefore improbable that the treatment of the cc binding should require large corrections. 
This conclusion is supported by the fact that the relative rate of ip' and direct J/ production 
(Eq. |l]), which at high energies should be independent of the production mechanism, is in 
agreement with experiment. 

In a leading twist description, an incorrect normalization of the charmonium production 
cross sections can arise from large higher order corrections or uncertainties in the parton 



distributions |]TT[. Taking into account that the normalization may be wrong by as much as 
a factor of 10 and that even such a /T— factor does not explain the polarization data of J/ip, 
a more likely explanation may be that there are important higher-twist contributions to the 



production of the J/ip and Xi as suggested in large xp case [EBl P7| 



Further theoretical work is needed to establish that the data on direct J/ijj and Xi 
production indeed can be described from higher twist mechanisms. Experimentally, it is 
important to check whether the J/ip^s produced indirectly via X2 decay are transversely 
polarized. This would show that X2 production is dominantly leading twist, as we have 
argued. Thus, the polarization of J/ip production from different channels provides a very 
sensitive discriminant of different production mechanisms. 

References 

[1] M. Vanttinen, P. Hoyer, S. J. Brodsky and W.-K. Tang, in preparation. 

[2] V. Barger and A. D. Martin, Phys. Rev. D31, 1051 (1985). 

[3] A. G. Clark, et al, Nucl. Phys. B142, 29 (1978). 

[4] R806: C. Kourkoumelis et al, Phys. Lett. B81, 405 (1979) 

[5] WAll: Y. Lemoigne, et al, Phys. Lett. B113, 509 (1982). 

[6] E673: S. R. Hahn, et al, Phys. Rev. D30, 671 (1984); D. A. Bauer, et al, Phys. Rev. 
Lett. 54, 753 (1985). 

[7] F. Binon, et al, Nucl. Phys. B239, 311 (1984). 

[8] E705: L. Antoniazzi, et al, Phys. Rev. Lett. 70, 383 (1993) 



7 



[9] E705: L. Antoniazzi, et al, Phys. Rev. D46, 4828 (1992) 

[10] E672: A. Zieminski, et al, Proc. XXVI Int. Conf. on High Energy Physics, Dallas, 
Texas, 1992, AIP Conf. Proc. No. 272, Ed. by J. R. Sanford, p. 1062. 

[11] G. A. Schuler, preprint CERN-TH.7170/94. 

[12] NA3: J. Badier, et al, Z. Phys. C20, 101 (1983). 

[13] C. Biino, et ai, Phys. Rev. Lett. 58, 2523 (1987). 

[14] E537: C. Akerlof, et al, Phys. Rev. D48, 5067 (1993). 

[15] E772: D. M. Aide, et ai, Phys. Rev. Lett. 66, 133 (1991). 

[16] J. F. Owens, Phys. Rev. D30, 943 (1984). 

[17] J. F. Owens, Phys. Lett. B266, 126 (1991). 

[18] R. Baler and R. Riickl, Z. Phys. C19, 251 (1983). 

[19] R. Gastmans and T. T. Wu, The Ubiquitous Photon: Helicity Method for QED and 
QCD, Clarendon Press, Oxford, 1990 

[20] J. Badier, et al, Z. Phys. CIS, 281 (1983). 

[21] J. S. Conway, et ai, Phys. Rev. D39, 92 (1989). 

[22] P. Chiappetta and M. Le Bellac, Z. Phys. C32, 521 (1986) 

[23] J. G. Heinrich, et al, Phys. Rev. D44, 1909 (1991). 

[24] W. Kwong, J. L. Rosner and C. Quigg, Ann. Rev. Nucl. Part. Sci. 37, 325 (1987). 
[25] L. Kopke and N. Wermes, Phys. Rep. 174, 67 (1989). 

[26] S. J. Brodsky, P. Hoyer, A. H. Mueller and W.-K. Tang, Nucl. Phys. B369, 519 (1992). 
[27] P. Hoyer, M. Vanttinen and U. Sukhatme, Phys. Lett. B246, 217 (1990). 



8 



This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9408371vl 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9408371vl 



